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^ ' Abstract 



The excess of high Q 2 events at HERA reported in the early 1997 by HI and 
ZEUS collaborations has become the subject of extensive studies in the framework 
of several models related to new physics. Here we concentrate on the most promis- 
ing, from our point of view, model describing HERA anomaly. We update our 
previous analysis || and take into account new HERA statistics of the 1997 year. 
HERA events are considered within the i?-parity broken SUSY model for a spe- 
cific scenario with several non-zero couplings. i?-parity broken SUSY with several 
non-zero couplings could explain both high Q 2 e + + jets and u + + jets observed at 
HERA. The consequence of such a particular scenario is the excess of high Q 2 di- 
or tri-jet events at HERA. The relation of this scenario for LEP and TEVATRON 
colliders is considered. This study shows that if a squark resonance does take place 
at HERA, supersymmetry with broken i?-parity can be revealed at either LEP200 
or TEVATRON in the near future. 
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1 Introduction 



Anomalous high-Q 2 events observed at HERA by HI JIJ and ZEUS |2| collaborations in 
the early 1997 have provoked much efforts of their explanation within different theoretical 
frameworks. For the 1994-96 year data an excess of high Q 2 events was at about 3a level. 
Later, these two collaborations reported an analysis with the 1997 year data included at 
the Symposium on Lepton-Photon Interactions |3[] and very recently at DESY seminar on 
the 13th of March, 1998 §. 

The present statistics is more than doubled (46.6 pb _1 - ZEUS and 37.04 pb _1 - HI) 
in comparison with the statistics of the 1994-96 years (20.1 pb _1 - ZEUS and 14.2 pb _1 - 
HI). The acceptance has also been increased. At present, for the 1994-97 data combined, 
the excess is at the 2a level, for example, for the above-mentioned luminosity the number 
of observed events with Q 2 > 20000 GeV 2 is 7 (ZEUS) + 10 (HI) = 17, while the number 
of predicted events is only 5.3 (ZEUS) + 4.4 (HI) = 9.7 

Basically there are three kinds of proposed explanations: 

1) contact-term interactions of very heavy bosons, heavy leptoquark-type objects being 
exchanged between a positron and a quark |6| ; 

2) non-standard evolution of the parton densities 0; 

3) leptoquark || or squark || resonance within the i?-parity violated SUSY model. 
The first two scenarios, as it has been shown, do not give quite a good explanation of 

HERA events jED]. 

We concentrate here on the possibility of the resonant e + q production at HERA which 
gives the most reasonable explanation of high Q 2 events. One can expect resonance 
production of LQ from valence ( u and d) or sea (u, d, s) quarks. 

However, the resonance production of e + u and e + d is excluded by HERA data for 
the e~p mode with 1 pb _1 statistics, because such signals would be seen in the e~u and 
e~d resonance production for quite a big coupling constant value of an order of 0.2 [fLU 
appearing in the scenario with a resonance in e + p collisions. 

On the other hand, the generally accepted leptoquark framework of Ref. [12] with 
a non-zero coupling only within one family is excluded for HERA by the combined 
CDF+D0 limit on the leptoquark mass of 240 GeV at 95% C.L. |L3|. This limit is 
established under the assumption that BR(LQ — > e + q) = 1. But for HERA, within 
this model, the branching ratio is fixed to unity if one takes into account the severe limit 
(5 x 10~ 7 ) on the product of chiral couplings from the branching ratio of the leptonic pion 
decay [[14 . 



That is why we would like to consider a more preferable, from our point of view 
scenario, when leptoquark appears as a scalar superpartner of a quark within the im- 



parity violating (RPV) SUSY model |15 |. 



We would also like to update our previous paper || with the 1997 HERA data included, 
namely, we consider a possibility of the i?-parity violation provided s-channel resonance 
at HERA for the scenario with several non-zero .R-parity violating couplings takes place. 

The starting point for our analysis is the superpotential of the Minimal Supersymmet- 
ric Standard Model with i?-violating terms included 

W = h E LE c H 1 + h D QD c H 1 + huQU c H 2 + fiH 1 H 2 

+f, i H 2 L i + \ ijk UL 3 El + \' ijk LiQjD c k + X'! jk U^D c k (1) 

Couplings Xijk and X"j k have the following properties: = —Xjik, A^- fe = — X" k j. 
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The terms with A and A' violate a lepton number L while the last term with A" violates 
a baryon number B. The B and L violating operators, if they exist simultaneously, lead 
to a fast proton decay ||21|| . Thus, only one kind of operators is allowed. HERA events 
can be explained in the framework of the .R-parity broken SUSY model only if A' 7^ 0. 
That is why we hereafter consider the case with A" = 0. 

The A and A' -terms result in the following Lagrangian in component fields: 



L Ijt v ,x = ^ijkftiLekRejL + e jL e kR u iL + e* kL (is iL ) c e jL - (i +-> j)] F + h.c, (2) 

e-ihdkRUjh ~ u jL d kR e iL - d* kL (e iL ) c u jL \ + h.c. (3) 
The strictest limits on different single couplings are presented in Table [I] (see, e.g. 0, 



Tjj and references therein). 
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Table 1: The strictest experimental bounds on .R-parity violating couplings for the squark 
mass rh = 100 GeV 

The bounds denoted by a dagger (t) are 2a bounds, the other bounds are at the la 
level. The bounds denoted by an asterisk (*) are based on the assumption of absolute 
mixing in the (SM) quark sector 
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The bounds in Table [l] are obtained from the following processes: 
[a) charged current universality [18, 19], 



BE 



(b) T(r — ► evv)/T{r — > [wv 

(c) bound on the mass of u e p0| , pT 

(d) neutrinoless double-beta decay 12 

(e) atomic parity violation 
(/) D° - D° mixing 0, |7 

(g) R e = T had (Z G )/T t (Z°) 

(h) r(vr -> eu)/T(n -> fiu) fl 
(1) BR{D+ -> K°*fi + ^)/BR(D + 
(j) v ii deep- inelastic scattering |T8| 
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(k) BR(T -> 7TZ/ r ) 1 

(I) heavy nucleon decay |J3 
(m) n — n oscillations |S2| 
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Bounds shown in Table |T| scale with the following factors: 
6), e), g) to to) - Af^/100 GeV; 

f) y/M-JIQO GeV- 



d) - (M§/100 GeV) 2 sjM- g /l TeV. 

Within the model under study only the following processes survive if one takes into 
account the above-mentioned low energy limits: 

e + d -> t, e + d -> c, e + s -> t (4) 

We would like to concentrate ourselves on the scenario with only a valent quark in- 
volved. The branching ratio of stop (scharm) decay into an electron and a quark should 
be less than unity in order not to contradict D0 and CDF data. There are two ways of 
having this: 

1) squark has also -R-parity conserving decay channels; 

2) squark has no -R-parity conserving decay channels but has additional R-parity violating 
decay channels, which leads in turn to the flavor changing neutral current (FCNC). The 
.R-parity conserving decay channels will be closed if, m 4 ~ < to& + and toj < raw + m~ b . 
The same situation could be for c-squark but for a smaller SUSY parameter space ||11|| . 

The first possibility has been intensively studied (see, for example, Ref. Jll], |35||). The 
.R-parity conserving decay will lead to the charged current event excess. Actually, HERA 
has some excess of such events [|J. Since the observed charged current events are events 
with a single jet and -R-parity conserving decays lead to multiparton final states, a special 
assumption is required that some partons should be too soft to be detected as a jet. 

The second possibility to have BR(q — > e + , d) < 1 has been considered in our previous 
paper ||. We considered the case of several non-zero couplings of the same order ~ 0.04. 
However, when one considers such a scenario, one should also take into account limits 
from the low energy experiments on the products of -R-parity violating couplings. Almost 
all combinations of several non-zero couplings contributing to the flavor changing neutral 
current are excluded. Modern limits on the products of -R-parity violating couplings come 
from: 

1) v — e conversion in nuclei ||36|| ; 

2) limits from rare leptonic decays of the long-lived neutral kaon, muon,tau and from 
mixing of neutral K- and B-mesons 



3) from B° decays to two charged leptons [[}§[ 

4) from muon decay to electron and photon ]39 

5) from features of K— and B— systems [10 . 

In Table |2] we present the most important for our analysis bounds from the literature 
cited above. 

After careful analysis of the relevant bounds only the following combinations of non- 
zero A' survive: 

1) for the c-squark resonance scenario: A' 121 , A 323 , A 333 , A 311 and A131; 

2) for the t-squark resonance scenario: A' 131 , A 333 , A' 311 and A131. In addition to these non- 
zero couplings there is a possibility for A' 112 or A' 113 to be non-zero as well. We will consider 
this case in relation to possible signals of -R-parity violation at LEP. 

We chose non-zero A"s equal to each other. An interesting fact is that bounds on the 
products of different A' close the decay channel both for t and c into a muon and a quark 
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and the only possibility for an additional channel of squark decay into r and b-quark is 
left (|A' 121 A' 222 |, |A' 131 A 232 | < 3.5 • 10" 7 , |37|). 

The purpose of the present paper is to consider this scenario with several non-zero 
i?-parity violating couplings for the explanation of the excess of high Q 2 neutral current 
DIS events at HERA as well as FCNC events observed by HI and discuss a possible check 
of such a scenario at LEP200 and TEVATRON. 

All analytic and most of the numerical calculations have been performed with the help 
of the CompHEP software package |5l[].This package allows one to make complete tree 
level calculations in the framework of any fed model. We have implemented into this 
package the part of the supersymmetric Standard Model with i?-parity violating terms 
relevant to our analysis. 



2 HERA events 

To describe HERA anomalous events, we used '94-97 data. For example, for Q 2 > 20000 
GeV 2 we have 17 events while only 9 are predicted by the SM. We have also taken into 
account the integrated luminosity of ZEUS+H1 equal to 83.6 pb -1 and an average detector 
efficiency 0.85. 

Distribution for Q 2 > Q 2 min for the Hl+ZEUS 1994-97 data and that for the SM 
prediction is shown are Fig. 0. One can see that the systematic deviation from the SM 
starts from Q 2 > 15000 GeV 2 , which is shown in Table [| 
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Figure 1: Q distribution, predicted by SM (solid line), SM+i?-parity breaking squark 
resonance (dashed line) [A' = 0.038] and data (solid dots). 



Combined systematic and statistical errors for a number of events, predicted by the 
SM are of an order of 10%. 

In Table ^ we present values for invariant mass of the events with Q 2 > 20000 GeV 2 . 



Table 4: 
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The events from the '97 year data are marked by an asterisk (*). We have taken mass 
values from the figure of 0. 

One can see that the events observed by the HI collaboration (Q 2 > 20000 GeV 2 ) 
are more concentrated around the invariant mass of 200 GeV than the ZEUS events. In 
conclusion of PH'97 it was said that ZEUS and HI data did not fit well together for the 
scenario with a single narrow resonance ||. But the statistics and mass resolution are 
not enough at the moment to exclude completely the scenario with one resonance: for 
Q 2 > 20000, 7 HI events and 3 ZEUS events are in the mass window 200 ± 20 GeV. 
Moreover, one can obtain a broad mass distribution if a superposition of two close mass 
eigenstates is considered 



In our calculations we have used the CTEQ3M structure function and Q scale 
equal to (p e ~ p' e ) 2 - The total cross section describing the observed data turned out to be 
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0.30 pb|, while the Standard Model gives only 0.16 pb. 

We use the set of non-zero i?-parity violating Yukawa couplings, which has been de- 
scribed above. We use: 

1) A 121 = A 323 = A 333 = A 311 = A', A131 = A 
or 

2) A 131 = A333 = A' 311 = A', A131 = A. 

We have calculated the complete set of Feynman diagrams for the process e + + p — > 
e + + jet + X in the framework of the i?-broken supersymmetric model for HERA energy. 
The relevant Feynman diagrams are shown in FigfJ. 



d d d d 




Figure 2: Diagrams for the e + + p — > e + + jet + X process in the supersymmetric model 
with i?-parity violation. 



As it is naturally expected, the main contribution comes from the diagrams with a 
resonant squark in the s-channel (99%). Contribution from the t-channel diagram is 
negligible (1%) for the case under consideration. However, these diagrams give the same 
order contribution as a nonresonant s-channel diagram with a heavy sparticle, but this is 
not the subject of our study. 

In our calculations we assumed NLO corrections by introducing the .fT-factor equal to 
1.25 for rrig = 200 GeV |5|. 



The presence of several non-zero couplings leads to the branching ratio BR(c(t) — > 
e + d) equal to 50%. 

We obtained the value of A' equal to 0.038 which describes the event excess in HERA 
data with the 200 GeV s-channel squark resonance. Fig. [I] shows the Q 2 distribution for 
the Standard Model and SUSY model with broken i?-parity tuned to describe data and 
the combined ZEUS+H1 data distribution itself. 

One can see that the Q 2 shape of deviation from SM is well described by the s-channel 
squark resonance. 

The presence of another non-zero coupling — A 323 or A 323 causes the flavour changing 
neutral current events: e + d — > rb. 

Leptonic and hadronic r-decay channels lead to the following signatures: 

1) e+ or fx + +If T +jet; 

2) di- or tri-jet events + Ife. 

Two outstanding FCNC events with a high pt fi + + jet(s) + missing Pt signature 
have been observed by HI collaboration and reported at Jerusalem conference ||42|| . Only 
0.05 events were predicted by the SM in the kinematic region of these two events. The 
scenario with non-zero couplings describes well the kinematics of the observed events. If 



3 We obtained this value using data presented by HI and ZEUS collaborations ||, g|. This value is 
the same as was presented at 'LP97: the statistics has been increased since that time and 3 new events 
for Q 2 > 20000 GeV 2 were observed by HI collaboration. 
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one takes into account the r — > ^ii7^v T decay branching ratio equal to 17%, then one will 
have about 2 fi + + jet + missing Pt events which were observed by HI collaboration. If 
the deviation from the SM for di- or tri-jet + If? events is found, it will be an additional 
favour for the scenario with several non-zero RPV couplings. 

3 Implication for LEP200 and TEVATRON 

Having in mind that processes with i?-parity violation could occur in e + e~ collisions at 
LEP, we consider the possibility of detecting it in the pair jet production via a sparticle 
exchange. Deviation of the jet production rate from the SM at LEP due to i?-parity 
violating SUSY has been the subject of several papers [||| §, (47j] . 

Comparing with the Standard Model, some excess in dijet production is expected. For 
this purpose, we have calculated the complete set of relevant diagrams presented in Fig.0. 
We consider the case with also non-zero A' 113 = A'. Hereafter, we assume that sleptons 
have a mass 200 GeV unless otherwise noted. We also assume that sleptons are not 
produced resonancely, considering as an example the scenario with -\/s e+e _ = 184 GeV. 

We would like to stress that there is a difference (which was missed in some earlier 
papers) in the up and down quark production at LEP for the scenario under study. For 
bottom-quark production, for example, there is an s-channel diagram with a sneutrino 
and t-channel diagram with an w-squark (Fig|| (a)). It makes the total cross section 
dependent not only on A' but also on A values. The t-channel diagram has negative 
interference while the s-channel diagram has a positive contribution to the total cross 
section. The total effect from both the diagrams is positive, when A' and A are of the 
same order of magnitude. At the same time, if only the t-channel diagram with .R-parity 
violating vertex is present, as it happens for the up-quark pair production (Fig.| (b)), 
the total effect is negative due to the negative t-channel diagram interference and the 
cross section is smaller than the SM one. This is illustrated in Fig.|4] for the 6-quark pair 
production (a) and the c-quark pair production (b), where the total cross section as a 
function of A' is shown (A is fixed equal to 0.038). 

If we take A131 equal to its upper limit (0.12), we have a 5% deviation (~ 0.21 pb) of 
the total cross section from the Standard Model for = 184 GeV LEP. In our calculation 
we put the value of A' to that describing the HERA data, namely, 0.038. If we assume 
that the b-jet production ratio can be measured with the 1% accuracy, one can establish 
the limit on A about 0.055 (A' = 0.038) or the limit on the product |A • A'| < 2 • 10~ 3 . 

As for the up-quark production, the t-channel diagrams bring a small negative contri- 
bution. One could expect some contribution of a 200 GeV SUSY particle to the Rb value 
very accurately measured at LEP90 which is equal to 0.2170 ± 0.0009 ||48[| . But it turns 
out that for yfs = 90 GeV the contribution from additional diagrams to the bb production 
is only in the 6-th digit. 

Though limits on A's are not so strict, it seems unnatural if their values are much higher 
than A'. We would like to look at the case when A is small or equal to and the excess in 
the bb production is absent. For example, for A ~ A' = 0.038 the bb production deviates 
from the Standard Model by less than 1%. By a lucky chance we have at present three 
different colliders with energy of an order of 1 TeV: lepton-lepton — LEP, lepton-hadron 
- HERA and hadron-hadron — TEVATRON. These machines are highly complementary 
to each other. We show below that possible failure to find the i?-parity broken SUSY for 
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Figure 3: Diagrams for the e + e — > qq process in the i?-parity violated SUSY model. 
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Figure 4: Total cross section of e + e — > qq versus A' for 6-quark (a) and c-quark (b) pair 
production at LEP. Dashed line — the SM prediction. 



our scenario (if it really exists) at LEP leads to the luck at TEVATRON in this search 
and vice versa. One of the most interesting effects for TEVATRON revealing the -R-parity 
violation is the single top production via a slepton decay |3|, [5(| which is resonancely 
produced in pp collisions. 

The Feynman diagrams for the 2^2 process of single top production are shown in 
Fig.| 

In comparison with our previous paper we take into account also 2^3 diagrams of 
single top production via i?-parity violating vertices: f-gluon fusion process and 2^3 
process with resonant f production. The corresponding diagrams are presented in Fig. [(]. 

In spite of the fact that the cross section is directly proportional only to A', there 
is also dependence of the cross section on the A coupling, which is related to the width 
and branching ratios of the selectron decay. For comparatively large A (0.1-0.2) f mainly 
decays leptonically and the single top production is suppressed. In this case, the deviation 
in the bb production could be measured at LEP. On the other hand, if A is small or even 
equal to 0, there is a significant contribution to the single top production via processes 
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Figure 5: Diagrams for 2 — > 2 single top production processes at TEVATRON in the 
i?-parity violated SUSY model. 
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Figure 6: r-gluon fusion and resonant r 2 — > 3 diagrams for single top production 
processes at TEVATRON in the .R-parity violated SUSY model. 



mentioned above. 

In our calculations we have fixed A' 311 and A 333 equal to 0.038 and studied the depen- 
dence of the single top quark cross section on the values of A131. We assume m to „=175 



GeV, use the CTEQ3M structure function [|3 and use the slepton mass as a QCD scale. 
We have calculated the following processes: 

1) ud — > tb, du — > ib; 

2) ug — > dp — > tfo-u, % — * t&<i, dp — > t&u. 

The process 2^3 has about 20% contribution to the total RPV cross section and thus 
should be taken into account. For example, for A 13 i=0, the RPV cross section of the 
2^2 process is 0.32 pb, while the RPV cross section of 2 — ► 3 is 0.078 pb. 

For our analysis we have used NLO cross sections for the SM single top production 




Figure 7: Total cross section of the pp — ► tb(q) single top production process versus A 
(a) and total cross section of the bb production at LEP versus A (b). Dashed line shows 
the SM prediction (l' 3U = X' 333 = 0.038) 
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using the results of Ref.|3(m t =175 GeV, ^=2 TeV): a{tb)=0.88 ± 0.05, a(Wg = tqb + 
tq) = 2.43 ± 0.4. NLO cross section is about 20% higher than the LO one (for details 
see [f^|). One could expect the same order corrections for the RPV diagrams of the single 
top quark production. We used this factor in our analysis. 

We can see that the contribution of the RPV processes can be as large as 20% of 
the total single top production rate. It is illustrated in fig [?] where the dependence of 
this deviation versus A (A' is fixed to 0.038) is shown. We include the results of the 
analysis for single top background processes into our study (see [^3j). A careful kinematic 
analysis allows one to suppress significantly the background and makes signal/background 
ratio equal to 1/2-1/3 |53[] which in turn gives the possibility to measure the single top 
production cross section with the accuracy 19% and 7% for TEVATRON Run2 and Run3, 
respectively |52| . At the same time, the top- and bottom-quark invariant mass distribution 
analysis will be more sensitive to the RPV single top physics than the total cross section 
measurement. 

We have found that if one takes into account the 10% accuracy of the single cross 
section measurement, the lower limit on A131 can be established: A131 > 0.07 (A' = 0.038). 
On the other hand, LEP data will allow one to put an upper limit on A131: A131 < 0.055 
(A' = 0.038). We can see that by lucky coincidence LEP and TEVATRON complement 
each other in search for supersymmetry with broken i?-parity. The existence of the RPV 
physics will be checked at LEP and TEVATRON which together with HERA could prove 
the validity, or restrict, or even exclude the scenario with several non-zero RPV Yukawa 
couplings. 

Conclusions 

The excess of high Q 2 HERA events can be reasonably explained within the Supersym- 
metric Standard Model with broken i?-parity. 

New 1997 year data have decreased the statistical significance up to la level. It is by 
no means conclusive and one should wait for new coming statistics. 

We have suggested some specific scenario for the .R-parity violating SUSY which has 
several non-zero couplings. This scenario could be a reasonable explanation of two out- 
standing FCNC muon events observed by HI collaboration. At the same time, one can 
expect an excess of di- or tri-jet events. We have checked that the Yukawa couplings (as 
well as their products) responsible for processes under study do not exceed their experi- 
mental bounds. Processes with .R-parity breaking can also take place in e + e~ collisions at 
LEP200 and TEVATRON, which are complementary to each other. If i?-parity violation 
really takes place, as has been shown above, it can be revealed at either TEVATRON or 
LEP in the near future. 
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